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Abstract. We use the single-histogram technique to study the critical behavior of the three-state Potts 
model on a (random) Voronoi-Delaunay lattice with size ranging from 250 to 8000 sites. We consider the 
effect of an exponential decay of the interactions with the distance, J(r) = Joexp(— ar), with a > 0, and 
observe that this system seems to have critical exponents 7 and 1/ which are different from the respective 
exponents of the three-state Potts model on a regular square lattice. However, the ratio 7/1/ remains 
essentially the same. We find numerical evidences (although not conclusive, due to the small range of 
system size) that the specific heat on this random system behaves as a power-law for a — and as a 
logarithmic divergence for a = 0.5 and a = 1.0 

PACS. 6 4.60. Fr - 5.10.Ln - 5.20.-y - 5.50.+q 



1 Introduction 

The randomness in the lattice of statistical spin models 
, has been studied in order to access the effect of impurities 
' and dilutions over their critical behavior. It was conjec- 
I tured by Harris that the sign of the critical exponent 
of the specific heat, a, determines whether the system is 
affected or not by such a randomness. For positive values 
of a, the impure system should have a critical behavior 
different from the one of the pure system. For negative 
[ values of a, on the other hand, the critical behavior of the 
. system should be the same for both cases. In the marginal 
■ case of a = 0, one should not be able to draw any conclu- 
, sion about changes in the system. 

The pure ferromagnetic three-state Potts model has 
I a = 1/3, hence, from the Harris criterium we expect to 
. find a different behavior with a random interaction sys- 
' tem. However, Picco used this model with two types of 
, interactions (Jo = 1 a-nd Ji — 1/10) randomly and equally 
distributed {p = 0.5) and did not find significant differ- 
ences from the pure case. In another study a quenched 
random interaction between nearest neighbors has been 
introduced, where the coupling factor Jy for each pair i,j 
is selected from two positive values J and J' with respec- 
tive probabilities equal to p and 1 — p. It was found that 
the exponent r]{= 2 — j/v) does not depend on the dis- 
order length, while and 7 vary continuously with this 
type of disorder, satisfying the concept of weak universal- 
ity proposed by Suzuki Q . The finite size scaling behavior 
of the Potts model with three and four states on regular 
square lattice has been studied by Kim and Landau S, 



where they found that for g = 4 the multiplicative loga- 
rithmic correction is insufficient to correct the dominant 
terms. 

In the present work, we investigate the effect of quenched 
disorder on the three-state Potts model using a two dimen- 
sional Voronoi-Delaunay network. This kind of disordered 
lattice exhibits a random coordination number that varies 
from 3 to 00, depending on the number of sites. In ad- 
dition, the distance r between nearest neighbors changes 
randomly from site to site. This geometrical feature is in- 
corporated in our model by using a coupling parameter 
that is an exponential function of r, J(r) oc exp(— or), 
with a > 0. In a previous work, we observed that the 
two-dimensional ferromagnetic Ising model with the same 
coupling mechanism and the same lattice displays criti- 
cal behaviors for a, (3 and v that are independent of the 
value of a and are the same as those of the pure system 

Here we extend this analysis to the three-state ferro- 
magnetic Potts model. 



2 Model and Simulation 



The Voronoi construction or tessellation for a given set of 
points is defined as follows. N points are randomly placed 
in a square with side lengths equal to N^^"^. Then, for each 
point we determine the polygonal cell consisting of the re- 
gion closer to this point than to any other. The points 
whose cells share an edge are considered neighbors. The 
dual lattice (Delaunay's) is composed by the vertices of 
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these polygons. The lattice obtained by linking the neigh- 
bor sites is the Voronoi network. 

The Hamiltonian for the three-state ferromagnetic Potts 
model is given by 



E 

<i,j> 



JljS{(7i — Uj), 



where K = l/ksT, T is the temperature, ks is the Boltz- 
mann constant and S is the Dirac's delta. The summation 
is performed over all pair of neighbors in the Voronoi's 
network and the spins can take values in the set {1, 2, 3}. 
We assume that 



Jin - Tj) = Joexp(-a|rj - r^l), 



where Jo is a dimensionless and positive interaction factor, 
ri and Vj are the position vectors of the sites i and j, 
respectively, and a > is a model parameter. 

For each grid size, we use the single-cluster algorithm 
llTll to simulate the system in the vicinity of the phase 
transition, i.e., with K — Kc-^q^, to determine the sys- 
tem finite size behavior. 

The value of ^Cmax determined in the following 
way. We find the first estimate Kq as the point of max- 
imum specific heat obtained in the curve (specific heat 
against K) constructed by simulations with K varying. 
The value of Kq is used to perform from B = 10 to 20 
independent simulations with 10^ MC steps each, from 
which, using the single histogram algorithm [ p^ , we con- 
struct a series of curves of C against K. We then average 
over all of these curves in order to find the specific heat 
curve as 



C(X) = (l/i?)^C«(if), 



(1) 



where C'^^^K) ^ K^N{< 

>K — < e >^) is computed 
using the single- histogram with the results of simulation i. 
Finally, we obtain, from the average curve, the maximum 
value of the specific heat, Cmax, and its abscissa ^Cmax- 
In order to estimate the ratios f3/v and j/iy, we use the 
fact that the magnetization at the inflection point and the 
susceptibility should scale, respectively, as 

I < Af > lij^f = L-^^'^fitL'/n cx L"'='/^ (2) 

Xmax{L) = x{Kc^^^iL),L)^AL'>^'', (3) 



and that 



Hence, the logarithmic derivatives should scale as 



dK 



In < lAf I > 



max 



_1 d_ 

\M\dK 



< \M\ > 



oc 



max 



and 
_d_ 

Ik 



In < |M1 > 



1 d 



|A'/2| dK 



(4) 
(5) 



max 



and so we can get a good estimate of ly from their scaling 
behavior. Finally, from finite size scaling arguments we 
can also predict that the specific heat should scale as 



Cmax(L) =Bo + BiliiL. 



3 Results and Conclusions 



(7) 



We study the critical 
three values of a (a = 

of a, we apply the finite size scaling technique |10 



behavior of the Potts model for 
0.0,0.5 and 1.0). For each value 

to- 



(6) 



gether with the single-histogram algorithm. We perform 
the same procedure for systems with different number of 
sites N = 250, 500, 1000, 2000, 4000, and 8000. The crit- 
ical temperature for infinite size system is estimated by 
using the fourth-order magnetization (Binder) cumulant 
and we find the critical values Kc = 0.607, 1.035 and 
1.959 and U* = 0.606, 0.615 and 0.623, corresponding to 
a = 0.0, 0.5 and 1.0, respectively. 

In Fig. 1 we show a log-log plot of M against L{= 
TV^/^) for a — 0.0, 0.5 and 1.0. By linear fitting each of 
these plots and using Eq. (||), we obtain P/v = —0.133, 
—0.118 and —0.106, respectively. The errors in these mea- 
surements are in order of ten percent, which is true also 
for the plots of the other figures. This is a consequence of 
the small range of the system size. 

In Fig. 2 we show the plot of xm against L also in 
logarithmic scale, from which we obtain the values of the 
exponent -f/v = 1.764, 1.751 and 1.754 for a = 0.0, 0.5 
and 1.0, respectively. As in Q, we observe that both 7 
and v vary with a while the ratio 7/1^ remains essentially 
the same. 

Figures 3 and 4 show the plots of the logarithmic deriva- 
tives (|^) and (H) against InL. The slopes of the curves pro- 
duce the estimates for l/v, from which we get the values 
1^1 = 0.840, 0.934 and 1.060 and 1^2 = 0.841, 0.934 and 
1.061 for a = 0.0, 0.5 and 1.0 respectively. 

In Fig. 5 we show the plots of Cmax versus InL for 
values of a — 0.0, 0.5 and 1.0. We observe that the curves 
for a = 0.5 and a = 1.0 can be well fitted by a straight 
line while the curve for a = 0.0 can not. The least-squares 
fits to data give the estimates Bq = —0.818 and 0.956 and 
Bi = 1.464 and 0.506 for a = 0.5 and 1.0 respectively. 
Figure 5 also contains the exponential fittings of Cmax 
versus L. From this figure, we can conjecture that, in the 
case of a = 0.0, the specific heat behaves like a power-law 
of L. In order to provide a quantitative support for this 
argument, we present in Table 1 the sum of the square 
errors obtained by the linear regressions (first row) and 
by the exponential fittings (second row) of Cmax x InL. 

Table 2 contains our numerical results in a condensed 
form. From the analysis of this data, we conclude that the 
randomness, which is introduced here through the geom- 
etry of the Voronoi-Delaunay lattice, changes the critical 
behavior of the system. This can be clearly observed by 
the change in the values of 7 and v. The value of /3 is not 
affected by the randomness. The analysis of the behavior 
of the exponent a is not conclusive since it is difficult, due 
to the limited number of data points, to distinguish the 
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trend (power-law or logarithmic) in the curve of Cmax ver- 
sus L. As already mentioned, the fittings shown in Fig. 5 
provide some indication that the specific heat is a power- 
law of L for a = 0.0, and varies as a linear function of InL 
for a = 0.5 and 1.0. This conjecture needs to be checked, 
however, for values of L > 8000. 
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Table 1. Total square residuals of linear and exponential re- 
gressions of Cmax X InL data 



a 


0.0 


0.5 


1.0 


Linear 
Exponential 


1.262949e-01 
6.369170e-03 


6.017198e-03 
1.474494e-02 


7.392830e-04 
3.294589e-03 
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Table 2. Theoretical and computed values of the critical ex- 
ponents 





13 


7 


V 






2d-Potts(g = 3) 


1/9 


13/9 


5/6 


2/15 


26/15 


a = 0.0 


0.112 


1.482 


0.840 


0.133 


1.764 


a = 0.5 


0.110 


1.635 


0.934 


0.118 


1.751 


a= 1.0 


0.112 


1.859 


1.060 


0.106 


1.754 



Figure Captions 
Fig. 1 - Logarithmic plots of magnetization (ln(M)) 
versus ln(L) for o = 0.0 (circles), 0.5 (squares) and 1.0 
(diamonds) . 

Fig. 2 - Logarithmic plots of susceptibility (ln(x)) ver- 
sus In(iy) for a = 0.0 (circles), 0.5 (squares) and 1.0 (dia- 
monds). 

Fig. 3 - Plots of In (I rf^ln < \M >|j^ax) ^^^^^^^ 
for a = 0.0 (circles), 0.5 (squares) and 1.0 (diamonds). 

Fig. 4 -Plots of In (I 3^ In < \A'P\ >|^^ax) versus ln(L) 
for a — 0.0 (circles), 0.5 (squares) and 1.0 (diamonds). 

Fig. 5 - Plots of Cmax versus ln(L) for a = (circles), 
0.5 (squares) and 1.0 (diamonds). Also shown are the lin- 
ear (solid lines) and exponential (dashed-lines) fittings to 
the data. 



